Quest for high T c in layered structures: the case of LiB 
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Using electronic structure calculation we study the superconducting properties of the theoretically- 
devised superconductor MSl-LiB (LiB). We calculate the electron-phonon coupling (A = 0.62) and 
the phonon frequency logarithmic average ((oj); 09 = 54.6 meV ) and show that the LiB critical 
temperature is in the range of 10-15 K, despite the frozen-phonon deformation potential being of 
the same order of MgB2. As a consequence, LiB captures some of the essential physics of MgB2 but 
(i) the electron-phonon coupling due to a states is smaller and (ii) the precious contribution of the 
n carriers to the critical temperature is lacking. We investigate the possible change in T c that can 
be induced by doping and pressure and find that these conditions cannot easily increase T c in LiB. 
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INTRODUCTION 

The quest for superconductivity in layered structures 
has become the focus of intense research since the dis- 
covery of superconductivity in MgB 2 (T c = 39 K) 
The layered structure of MgB2 generates one of its most 
prominent features, namely the B 2p xy orbitals form a- 
bands 0, 0, Q which are weakly dispersing along the 
k z direction and have a marked two dimensional charac- 
ter. In MgB2 the a— bands are hole doped, but the top 
of these bands is only w 0.5 eV higher than the Fermi 
level. The a— bands Fermi surface sheets [3], two slightly 
warped cylinders with axis perpendicular to the boron 
layers, generate a huge electron-phonon coupling along 
the k z direction. The carriers in the 7r— bands, formed by 
the B 2p z orbitals, further enhance the average electron- 
phonon coupling [5(. 

The formation of a and ir states is typical of graphite- 
like structures composed by boron or carbon atoms. 
Given the success of MgB2 it is natural to look for high T c 
superconductivity in structures having similar features. 
The problem is that, given the boron layers, small vari- 
ations in the valence or mass of the intercalant or in the 
structural parameters are sufficient to considerably alter 
the a or tt bands positions or the shape of their Fermi 
surfaces and consequently destroy superconductivity. For 
one or some of these reasons AIB2, ZrB2, NbB2, M0B2, 
YB 2 , TaB 2 , TiB 2 , HfB 2 , VB 2 and CrB 2 are not super- 
conducting @, 0, Q • 

The hope of finding new superconducting materials in 
layered structures was recently increased by the discovery 
of superconductivity in the graphite intercalated com- 
pounds, YbC6 and CaC6 0, [H3]. This is particularly 
promising since a huge number of intercalants are avail- 
able for graphite ll|. In CaC6, despite the layered struc- 



ture and the existence of a and tt bands originated from 
the carbon 2p orbitals, the electronic structure close to 
the Fermi level is completely different from that of MgB2. 
The 7T bands, reminiscent of the graphite ones, and an in- 
tercalant free-electron-like band [12|, [13J cross the Fermi 
energy. The intercalant band forms a spherical Fermi 
surface [THI l . The electron-phonon coupling of CaC6 
(A = 0.83) is mainly due to coupling of the interlayer 
band with C vibrations perpendicular to the graphite 
layers and with Ca vibrations. So, even though miss- 
ing the c-bands, CaCg reaches an interesting 11.5 K T c . 
This temperature is substantially enhanced by pressure 
(T c = 15.1 K at « 8 GPa |l6[), contrary to what happens 
in MgB 2 . 

As can be seen from the above examples, even if one re- 
stricts to sandwich structures formed by boron or carbon 
layers, the details of the electronic and phonon spectra 
and, subsequently, the critical temperature can change 
dramatically when the intercalant is included. As a con- 
sequence, a theoretical approach is absolutely necessary 
to identify the most probable superconductors or at least 
to exclude the less probable ones. 

An attempt in this direction has been recently made 
in Ref. 17|, where by using ab initio methods the au- 
thors studied the possible hole-doping of LiBC, awl 
eV gap semiconductor. The authors suggested that a T c 
of the order of MgB2 could be reached if the insulating 
LiBC is substantially doped with holes to obtain Lio.sBC. 
Successive experimental studies have indicated that the 
structural response to the introduction of holes unfavor- 
ably modifies the electronic structure of Li^BC, and so 
far no high T c superconductivity has been found in this 
system [18j. 

Ideally designing new superconductors ab initio re- 
quires three steps. The first is the determination of the 
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most stable structures given a set of atomic species. The 
second is the calculation of the electronic structure to 
verify that the given structure is at least metallic or can 
be made metallic easily. The third is the determination 
of the phonon dispersion and of the electron-phonon pa- 
rameters. 

The first point is a daunting task even if one restricts 
one's search to a specific set of likely candidates [l9j. A 
systematic approach to tackle this problem has been re- 
cently offered in the way of data mining of ab initio cal- 



ture basing our choice on the following considerations. 
On the one hand, it has the smallest unit cell of all MS 
structures, which offers computational efficiency. On the 
other hand, even though other stacking sequences are 
possible (e.g. MS2, Refs. [H, [H), MSI is a good rep- 



culations [2fj, |21|, |22j. In this method one uses the in- 
formations obtained from ab initio calculations of many 
different structures to build a database that can be then 
used to judge the stability of new structures. Application 
of this method to intermetallics has led to the identifi- 
cation of new layered lithium monoboride phases which 
have a goo d chance to form under proper synthesis con- 
ditions [231 ]. 

Once a stable metallic structure is given, a calcula- 
tion of the phonon spectra and of the electron-phonon 
coupling needs to be performed to obtain T c . Indeed, 
while some qualitative information can be extracted from 
electronic-structure [2j, for a quantitative analysis step 
three is absolutely necessary. 

In this work we investigate the superconducting prop- 
erties of the previously determined metal sandwich (MS) 
lithium monoboride [23j by calculating its phonon spec- 
trum and electron-phonon parameters. This system is 
metallic and, from qualitative arguments, one can infer 
that T c is of the same order of that of MgB2 [5^. In- 
deed this system has an electronic structure which is a 
hybrid between those of MgB2 and CaC6, since there 
are hole-doped cr— bands forming cylindrical Fermi sur- 
faces and there is an intercalant band crossing the Fermi 
level. Moreover the deformation potential is compara- 
ble to that of MgB2 [5^1 . From this point of view, LiB 
is a much more promising material than Lio.sBC, be- 
cause even without doping it has a significant density of 
cr-states at the Fermi level. 

Unfortunately, the full electron-phonon coupling calcu- 
lations performed in this paper indicate that LiB should 
have a T c in a 10-15 K range. We show that LiB captures 
some of the important physics of MgB2, namely the role 
of the a— bands, but it lacks the contribution of the it 
states to the electron-phonon coupling and it is only a 
far relative of CaC6 because the interlayer band is very 
weakly coupled with the phonons. In an attempt to im- 
prove the situation we examine what role the hydrostatic 
pressure and doping can play in determining the critical 
temperature. 



TECHNICAL DETAILS 

In all our calculations of the layered lithium mono- 
boride we have used the MSI theoretical crystal struc- 



resentative model of the layered lithium monoboride be- 
cause the long-period shifts are expected to have little 
effect on its superconducting properties [25]]. MSl-LiB 
has a rhombohedral unit cell with R3m space group. 
There are four atoms in the primitive unit cell with 
Wyckoff positions Li(2c)(l/2 - z u , 1/2 - z u , 1/2 - z u ), 
Li(2c)(l/6 + z u , 1/6 + z Li , 1/6 + z u ), B(2c)(-S, -6, -S) 
and B (2c) (2/3 - d, 2/3 - 5, 2/3 - 5). The fully relaxed 
parameters arc a = b = c = 5.92 A, a = (3 = 7 = 29.8°. 

Density functional theory (DFT) calculations are per- 
formed using the Quantum Espresso code [26|_within 
the generalized gradient approximation (GGA) [27] . We 
use norm-conserving pseudopotentials [28| with configu- 
ration 2s 1 2p° and non-linear core correction [29] ] for Li, 
and configuration 2s 2 2p 1 for B. The wavefunctions are 
expanded using a 50 Ry cutoff. The dynamical matrices 
and the electron-phonon coupling are calculated using 
density functional perturbation theory in the linear re- 
sponse [26] . For the electronic integration in the phonon 
calculation we use an Nk = 12x12x12 uniform It- 
point mesh and Hermite-Gaussian smearing from 0.05 
Ry. For the evaluation of the electron-phonon coupling 
we use an Nk = 40 x 40 x 40 Monkhorst-Pack mesh. For 
the A average over the phonon momentum q we use an 
N q = 4 x 4 x 4 q— point mesh. The phonon dispersion is 
obtained by Fourier interpolation of the dynamical ma- 
trices computed on the N q mesh. 

The pressure- and doping-induced changes in the elec- 
tronic properties of LiB are studied with Vienna Ab initio 
Simulation Package VASP 0,[3l| within the GGA [27 1. 
We use projector augmented waves (PAW) [32| pseudopo- 
tentials, in which Li semi-core states are treated as va- 
lence; the energy cutoff is set at 30 Ryd. The projected 
electronic density of states (EDOS) is found by decompo- 
sition of the wavefunction within a sphere of the default 
PAW radius of 1.7 a.u. For the MSI unit cell, the 2 x 2 x 3- 
MS1 and 2 x 2 x 1-MS2 supercells we use 31 x 31 x 31, 
18 x 18 x 6, and 18 x 18 x 10 Monkhorst-Pack k-meshes, 
respectively. 



BAND-STRUCTURE, DOS AND FERMI 
SURFACE 

The band structure of Li 2 B 2 is presented in Fig. [1] 
(see footnote [H} for high-symmetry points notation). 
Similarly to what happens in MgB2 [2|, there are two 
boron <r-bands crossing the Fermi energy tp. Compared 
to the tr-bands in MgB2, these bands are even more two- 
dimensional (due to the larger interlayer distance) and 
shifted by more than 0.6 eV to higher energies at the T 
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FIG. 2: (Color online) Fermi surface of MSl-LiB. For conve- 
nience, the Brillouin zone is stretched along the z-direction. 



FIG. 1: (Color online) Band structure of LiB. The size of the 
empty (full) dots represent the amount of Li (B) character at 
a given k-point. See footnote [33|] for high-symmetry points 
notation. 



point. As in MgB2, they generate two cylindrical Fermi 
surfaces (in our case with axes along the FT direction, 
Fig. [2]). The boron ir states in LiB resemble more the 
7T states of graphite, as they cross exactly at ep, so that 
LiB is lacking 7r Fermi surfaces altogether. In MgB2 these 
states cross at about 2 eV above ep, which leads to the 
appearance of an extended 7r Fermi surface 0] ■ Another 
important difference between the electronic structures of 
the two borides is the presence of a lithium band at ep in 
LiB. The position of this band resembles the intercalant 
band in CaC6 although in LiB it has substantial 

hybridization to boron states close to the T-point. The 
corresponding Fermi surface (a compressed sphere) is de- 
picted in Fig. [2j 

The total density of states (EDOS) and the EDOS pro- 
jected over atomic orbitals is illustrated in Fig [3] The 
main component at ep is given by boron pa states. As in 
graphite the boron p7r EDOS at ep is zero and increases 
slowly and linearly immediately after ep. 



PHONON SPECTRUM AND 
SUPERCONDUCTING PROPERTIES 

The phonon dispersion and the phonon density of 
states (PHDOS) are illustrated in Fig. 0] The phonon 
modes at the T-point are decomposed as 1Ai g + lA^u + 
2E g + 2E U [34\. To distinguish between modes with the 




FIG. 3: (Color online) Electronic density of states projected 
over selected atomic orbitals. 



same symmetry but different eigenvectors we use the fol- 



lowing notation: Ai g [B 2 ], E g [B a 



K [(Li - 
[(Li + B) 



B) K y], 

,1, E+ 



A\ g [Li,], A- 2u [(Li-B),], E* g [LiJjAl 
[(Li + H) xy \, where in brackets we give the corresponding 
atoms and vibrations. For convenience, wc label phonon 
branches everywhere in the Brillouin zone using the name 
of their representation at F. 

Except for acoustic modes, a clear separation exists be- 
tween optical Li and B vibrations. Li modes are confined 
in the 40-55 meV region and are not dispersive, meaning 
that Li-vibrations behave essentially as Einstein modes. 
Boron in-plane vibrations are softened along the TT di- 
rection due to coupling to the a bands. The softening at 
r of the Eg phonon branches is approximately 20 meV, 



4 




♦ 


Lixy 




Li z 


□ 


Bxy 




Bz 



i 1 

f 



FIG. 4: (Color online) Phonon dispersion of MSl-LiB with 
decomposition in in plane Li and B vibrations (labeled Li xy 
and B xy respectively) and out-of-plane Li and B vibrations 
(labeled Li z and B z respectively). Notation for the phonon 
modes at the F-point, given next to the vertical axis, is ex- 
plained in the text. 



to be compared with the almost 40 meV in MgB2 for the 
E 2g modes [1, HH . This suggests a strong coupling of the 
a bands to the in-plane vibrations in LiB 23[ , but not as 
strong as in the case of MgB2. 

The three acoustic modes along the FT direction are 
substantially softer with respect to the other directions. 
At the zone border, T, these modes are formed by (i) in- 
plane (Li+B)^ vibrations at energies ss 8.5 meV and (ii) 
out-of-plane (Li+B) z vibrations at ss 8.0 meV. They can 
be related to the soft modes at the T point in MS2-LiB, 
discussed in Ref. [23[ . 

The superconducting properties can be understood cal- 
culating the electron-phonon coupling A qt/ for a phonon 
mode v with momentum q: 



A ql / — 



uj^N(0)N k 



\9kn, 



k+qm I 



2 £(ekn)£(ek+qm) 



k,n.m 



(i) 

where the sum is over the Brillouin Zone. The matrix 
element is g£„,k+ q m = (kn|5V/5u qw |k + qm)/^/2o; q)y , 
where it q „ is the amplitude of the displacement of 
the phonon, V is the Kohn-Sham potential and N(0) 
is the electronic density of states at the Fermi level. 
The calculated average electron-phonon coupling is A = 
Sqiy ^qy/N q ~ 0.62 (JVjfe and N q are the previously 
defined k-space and q-space mesh dimensions, respec- 
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FIG. 5: (Color online) Phonon density of states (PHDOS), 
partial phonon density of states projected over selected vibra- 
tions, Eliashberg function and integrated Eliashberg function 
for MSl-LiB. For clarity, the integrated Eliashberg function 
has been multiplied by a factor of 4 



tively). 

The Eliashberg function 



a 2 F{uj) = — Y A q „u ql ,<S(w - w q „) 



(2) 



and the integral A(cj) = 2 du' a 2 F (lu 1 ) / oj' are shown 
in Fig. [5J As can be seen most of the contribution comes 
from phonon states in the 60-90 meV region and a smaller 
contribution comes from low energy states. 

An estimate of the different contributions of the in- 
plane (Li xy and B^) and out-of-plane (Li z and B 2 ) vi- 
brations to A can be obtained from the relation 



\jf3,ia 



(3) 



iaj/3 



iaj/3 



where i, a indices indicate the displacement of the 
i th atom in the Cartesian direction a, [Gq]i a ,j/3 — 
J2k, n , m A 9* a 9j0S(^n)S(e k+ q m )/[N(O)N k }, and g ia = 

(kn\5V/Sxqia\k + qm)/\/2. The C q matrix is the 
Fourier transform of the force constant matrix (the 
derivative of the forces with respect to the atomic dis- 
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placements). The decomposition leads to: 
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(4) 



The off-diagonal terms are small (but not negligible) 
compared to the total A. Most of the coupling is to the 
in-plane B vibration; contributions from the Li and the 
out-of-plane B vibrations are smaller. Since the cr-bands 
do not couple to the B z vibrations and since there are no 
7r Fermi surfaces, the coupling to B z vibrations is due to 
the intercalant band. Note that the decomposed values 
of A contain contributions from different modes and are 
summed over all the q-points in the Brillouin zone. For 
example, Ab xb ,b xh = 0.46 includes the coupling to the 
in-plane E g , E~, and E£ branches. By examining the 
integrated Eliashberg function X(uj) in Fig. 03c) and the 
phonon characters in Fig. [4] one can infer that the E g 
branch is the most important of the three: among them 
it has the highest PHDOS in the 70-100 meV range, in 
which A gains most of its total value. The soft in-plane 
E+ branch is far less important, as the net contribution 
from all the soft modes having energy under 20 meV is 
only « 0.08 (Fig. 02(c)). 

It is instructive to compare our result with other lay- 
ered superconductors. In MgB2 the coupling of the cr- 
bands to the phonon modes is X^ ! § B ' 2 = 0.62 ± 0.05 
[Bj], while in LiB the corresponding value is less than 
0.46, as discussed above. This difference can be clari- 
fied by noting that the E 2g phonon linewidth 7 qi _e 2 — 
27r./V(0)a; q £ 2 A ql/ along TA in MgB2 happens to be com- 
parable in magnitude with that of the E g mode along FT 
in LiB. Therefore, the reduced electron-phonon coupling 
in LiB is mainly due to the E g phonon frequency being 
harder than the E 2g one in MgB 2 (to^f /uj^f 2 « 1.3 
at L). This unfortunate result can be linked to the ab- 
sence of the 7r carriers, which play an important role in 
softening of the E 2g mode in MgB 2 3(| 37 1. 

We find that LiB and graphite intercalated compounds 
have few similarities in terms of superconducting fea- 
tures. In particular, in CaCg the intercalant modes 
are responsible of ~ 50% of the total electron-phonon 
coupling, and the rest comes from vibrations of carbon 
modes in the direction perpendicular to the graphite lay- 
ers. In CaC6 one has Aca xy + Aca z = 0.33 and Ac z = 0.33 
[lH . In LiB the overall contribution of B z , hi xy and Li z 
vibrations is less than half of that of CaCg, which means 
that while LiB captures some of the physics of MgB2, 
it does not capture the physics of graphite intercalated 
compounds to full extent. This is also clear from the 
phonon spectrum of CaCe where the intercalant modes 
are at energies lower than 20 meV and one of the Ca 
modes undergoes a marked softening with a correspond- 
ing large electron-phonon coupling (at point X of Fig. 2 



in Ref. [13[). In LiB, on the contrary, the Li modes are 
much higher in energies (~ 50 meV) meV and disper- 
sionless. The main reason for this difference comes from 
the mass of Li which is 5.77 times smaller than that of 
Ca leading to frequencies which are on average 2.4 times 
larger. 

The critical superconducting temperature is estimated 
using the McMillan formula [38| : 



1.2 



exp 



1.04(1 + A) 



(5) 



A-/i*(l + 0.62A) 
where fi* is the screened Coulomb pseudopotential and 

(C0)l g = e !/0 +O °« 2 -FV)l°gM/-^ (6) 

the phonon frequencies logarithmic average. We obtain 
(^)iog = 54 meV leading to T c of approximately 10-15 K 
for /i* = 0.14 — 0.1. This value could be further enhanced 
by multiband effects. 

DOPING AND PRESSURE EFFECTS 

Even though a theoretically-devised from scratch su- 
perconductor with T c = 10 — 15 K could be considered 
a success of the materials prediction methodology, the 
stoichiometric LiB compound falls short of expectations 
to compete with the record- holding binary MgB 2 . In 
this section we investigate whether it is possible to favor- 
ably modify the electronic properties of LiB and achieve 
higher T c by doping or applying pressure. We pay spe- 
cial attention to the evolution of the tt states, since their 
reintroduction at ep may soften the E g mode and lead 
to a larger coupling. 

As has been pointed out previously [23|, the bonding 
7r states are completely filled under ambient conditions. 
Because the band crossing in LiB at the Fermi level is 
accidental, it may be possible to move the crossing point 
with pressure and increase the 7r-bands EDOS at 6f- Fig- 
ure 02a) reveals that there is indeed a rapid change in the 
a and ir EDOS, followed by a plateau after 5 GPa. This 
behavior is a reflection of two distinctly different regimes 
of the LiB structural changes: i) in the 0-5-GPa pressure 
range the Li-Li interlayer spacing quickly shrinks and the 
B-B bond length slightly expands so that at 5 GPa they 
become about 0.5 and 1.02 of their zero-pressure values, 
respectively; ii) for pressures above 5 GPa the Li atoms 
in the bilayer start experiencing the hard-core repulsion 
and the compound compresses more isotropically. The 
inset in Fig. [6fa) illustrates that by reducing the Li-Li 
interlayer distance one forces the charge from the inter- 
calant band (completely emptied at about 6 GPa) into 
the boron n and a states (lowered by 0.7 eV at that pres- 
sure). Once the charge redistribution is complete, no ap- 
preciable changes in the EDOS are seen for the boron 
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FIG. 6: (Color online) a) EDOS at the Fermi level projected 
over a and n states as a function of pressure; b) the same as a 
function of doping level: the hollow points represent charged 
LiB unit cell in a neutralizing background; the solid points 
correspond to substitutional Mg and Al doping of LiB in the 
MSI and MS2 superstructures described in the text [4lj]. For 
comparison, the a and tt EDOS in MgB2 are 0.098 and 0.064 
states/(eV-spin-f.u.), respectively. The inset shows pressure- 
induced changes in the position of different states near the 
Fermi level: cr-boron and intercalant (labeled £) states at the 
r point and crossing of the tt-tt* states along the TS direction. 



states up to at least 30 GPa. Therefore, the peculiar be- 
havior of the nearly free electron intercalant states (also 
observed in other systems [l2, 39]) is the only meaningful 
factor allowing modification of the LiB boron states with 
pressure. 

These simulations demonstrate that the compression of 
LiB does not lead to the desired 7r-bands EDOS values 
comparable to those in MgB2. Moreover, the hydrostatic 
pressure causes such a quick drop in the tr-bands EDOS 
that this will likely negate any possible enhancement in 
the electron-phonon coupling due to the reintroduction of 
the Tr-states at €f ■ The phonon modes are also expected 
to harden under pressure, further reducing the electron- 
phonon coupling in LiB [40j. 

LiB has plenty of available bonding a states, there- 
fore the compound should be easy to electron-dope. A 
quick examination of the boron EDOS states around the 
Fermi level (Fig. [3] or Fig. 4 in Ref. 23]) gives an 
idea on what possible changes in the Fermi surfaces and, 
eventually, in the electron-phonon coupling the doping 
could lead to. At small doping levels, when the rigid 
band approximation normally holds, the EDOS from the 
two dimensional a bands should only slightly fluctuate 
until the states are completely filled, which happens at 
Aq « eN(0)A£ w 0.70 (e/eV) (1 eV) = 0.70 e/f.u. = 
0.35 e/boron. The EDOS from the tt bands grows slowly 
and even at the relatively high electron-doping of 0.35 
e/boron it would amount only to about a half of what 



is observed in MgB2. Another way to tweak the elec- 
tronic structure could be to hole-dope LiB as it is done 
for Li^BC 18]. The known limitations of this approach 
are the buckling of the hexagonal layers and the eventual 
dcstabilization of the compound upon heavy Li depletion 

We first simulate the electron- (hole-) doping using a 
charged cell with a neutralizing positive (negative) back- 
ground. Normally, in this approach one can safely relax 
the unit cell parameters and obtain valuable informa- 
tion about the bond length variation under small dop- 
ing. However, in the case of the electron-doped LiB the 
repulsion between the negatively charged boron layers 
overcomes the weak binding between the lithium layers, 
causing the c-axis to undergo unphysical expansion even 
at small levels of doping. Therefore, we fix the c-axis at 
the zero-doping value and relax only the remaining three 
parameters. The set of data, shown as hollow points in 
Fig. EHb), supports our earlier conclusion that the tt- 
band EDOS cannot be easily increased. Note that the 
approximations used in this test, i.e. the fixed c-axis and 
the use of a neutralizing background, may influence the 
results to some extent. For example, the positive elec- 
trostatic potential from ionized dopants could bring the 
delocalized Tr-states down (in addition to the rigid band 
downshift) and could potentially be an important factor 
in increasing the tt EDOS. 

To address these limitations we use a more realistic 
model of the electron-doped LiB by substituting Li with 
Mg or Al. Small doping levels are obtained only for large 
unit cells; we use the hexagonal 2 x 2 x 3-MS1 and 2x2x1- 
MS2 supercells with 48 and 32 atoms, respectively. Re- 
placement of one or two Li atoms in these structures re- 
sults in the 1/24, 1/16, 1/12, and 1/8 concentrations of 
dopants per boron, and the level of doping is found in the 
assumption that they give up all their valence charge. In 
all the cases the c-axis expansion in the fully relaxed unit 
cells does not exceed 6%. The resulting averaged boron 
EDOS for the Mg and Al sets are shown in Fig. [Sfb) 



as solid points [4JJ. The scattered presence of dopants 



in the lattice should cause some dispersion of the local 
boron properties. A general trend observed in our su- 
percell calculations is that a downshift of the tt and a 
states happens only for B layers in direct contact with 
the dopant. Typical values of the downshift that a single 
Mg (Al) atom induces in all eight atoms in a neighboring 
B layer are about 0.2 (0.5) eV. It is not easy to isolate the 
importance of different factors defining the level of B dop- 
ing, i.e. the simple charge transfer, the electrostatic effect 
discussed above and the structural changes (expansion of 
the c-axis and contraction of the B-B bond). However, 
Fig. [Hb) demonstrates that the net effect of the substi- 
tutional doping is described reasonably well within the 
rigid band model. 

In summary, our tests indicate that it is rather difficult 
to reintroduce a significant amount of tt states at cf with 
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hydrostatic pressure or small doping, because the band 
crossing in LiB happens to be exactly at cf, about 2 eV 
lower than in MgB 2 . To have a chance of substantially 
increasing T c , one should search for more radical ways 
of modifying the electronic structure of the MS metal 
borides. 



CONCLUSIONS AND PERSPECTIVES 

In this work we have investigated electron and phonon 
properties of the recently theoretically-devised supercon- 
ductor LiB [23| . By studying in details the phonon prop- 
erties of this hypothetical material we have found that 
its critical temperature is of the order of 10-15 K. Super- 
conductivity is mainly of the MgB2 kind with planar pa 
bands strongly coupled with phonons. Differently from 
MgB2, LiB has no Fermi surface generated by tt states 
but an additional intercalant one. Thus, its electronic 
structure can be seen as a hybrid between MgB 2 and 
CaC6- However, the intercalant electronic states of LiB 
are weakly coupled with the Li and B 2 vibrations so that 
the overall electron-phonon coupling is only A = 0.62. 

Since the discovery of MgB 2 , no other diborides have 
been found with high T c (for a full list, see Ref. [12] )• 
If we compare LiB with the known diborides, our calcu- 
lated 10-15 K critical temperature is not so low, although 
it is far from the 39 K of MgB2. Nevertheless, the study 
of LiB gives an important understanding. The common 
belief is that the main effect for the singular and unique 
behavior of MgB2 is given by the presence of almost two- 
dimensional a bands. LiB has even more planar a bands 
and even higher a EDOS at relative to that in MgB2 
23L 24 1; however, their contribution to the total electron- 
phonon coupling turns out to be at least 25% smaller. 
This reduction can be attributed to the differences be- 
tween the in-plane boron vibrations in the two borides, 
caused mainly by the lack of the n carriers at ep in LiB. 
Namely, the softening of the E g mode in LiB is substan- 
tially smaller (by about a factor of two) than that of E^ g 
in MgB2, which makes the former mode be noticeably 



harder (u;^ /w 



MgB 2 



1.3 at r). Therefore, in addition 



to the direct loss of the tt states contribution to A, their 
absence at cf in LiB also has a strong indirect negative 
effect on the overall electron-phonon coupling. 

We have investigated whether this somewhat unex- 
pected obstacle, preventing LiB to be a truly high T c 
superconductor, could be overcome with moderate mod- 
ifications of the compound's properties. Behavior of elec- 
tronic features important for the LiB superconductivity 
has been examined as a function of small doping and pres- 
sure. Our results indicate that the tt EDOS increases very 
slowly and cannot reach the desired values, at least not 
before the a EDOS is substantially reduced. Thus, small 
doping and pressure are not expected to significantly im- 
prove T c in LiB. A promising direction to fix the problem 



would be to find a suitable LiB-based ternary alloy; this 
question is currently under investigation. 

Theoretical development of potentially important su- 
perconducting materials is a difficult task because T c 
critically depends on their band structure features and 
vibrational properties. The challenge is even greater if 
one attempts to design a superconductor from scratch, 
since one first needs to ensure its thermodynamic sta- 
bility. The case of LiB shows that it is possible to the- 
oretically predict a compound that both i) has a good 
chance to form and ii) possesses interesting supercon- 
ducting properties. Study of such promising candidates 
gives important insights into how to perform a more tar- 
geted search for novel superconducting materials. 

While we were finishing writing this paper, a preprint 
on the related structure MS2-LiB appeared on-line 37 1. 
The results of the paper are similar to ours except for 
some numerical details that can probably be related to 
the different unit cells and k,q-point samplings usedQ. 
However, our conclusions concerning the possibility of 
increasing the T c in LiB by doping are rather different, 
as explained in the previous Section. 
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